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Abstract 



A sample of 7253 — 7r^e+e~(7) decay candidates with 1.0% background 
contamination has been collected by the NA48/2 experiment at the CERN SPS, 
which allowed a precise measurement of the decay properties. The branching ratio 
in the full kinematic range was measured to be BR = (3.11 ± 0.12) x lO"'^, where 
the uncertainty includes also the model dependence. The shape of the form factor 
W{z), where z = (Mee/Mx)'^, was parameterized according to several models, and, 
in particular, the slope 5 of the linear form factor = Wo{l+6z) was determined 

to be = 2.32 ± 0.18. A possible CP violating asymmetry of and K~ decay 
widths was investigated, and a conservative upper limit of 2.1 x 10~^ at 90% CL 
was established. 
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Introduction 



Radiative nonleptonic kaon decays represent a source of information on the structure 
of the weak interactions at low energies. The flavour- changing neutral current process 
TT^e+e^, induced at one-loop level in the Standard Model and highly suppressed 
by the GIM mechanism [T], is of particular interest. The K"^ 7r^/+/^ processes have 
been described by the Chiral Perturbation Theory (ChPT) [2]; several models predicting 
the form factor characterizing the dilepton invariant mass spectrum, and thus the decay 
rate, have been proposed [H HI [5]. 

The first observation of the — > vr+e^e^ process was made at the CERN PS more 
than 30 years ago [6] , followed more recently by BNL E777 [7] and E865 [8] measurements. 
The most precise of these, E865, based on a sample of 10300 candidates with 1.2% back- 
ground, allowed a detailed analysis of the decay form factor and rate, and a test of the 
next-to-leading order ChPT calculation |3]. 

In this paper, a new precise measurement of the —>■ n^e~^e~ decay is reported, 
based on the full data set collected in 2003-2004 by the NA48/2 experiment at the CERN 
SPS. In addition to spectrum and rate studies, the first observation of the K~ decay 
allowed setting an upper limit on the charge asymmetry of — > n~^e^e~ and 
7i~e~^e~ decay rates, which can be related to CP violation. 

1 The NA48/2 experiment 

The NA48/2 experiment, specifically designed for charge asymmetry measurements [9], 
uses simultaneous and K~ beams produced by 400 GeV/c primary SPS protons 
impinging on a beryllium target. Charged particles with momentum (60 ± 3) GeV/c 
are selected by an achromatic system of four dipole magnets with zero total deflection 
('achromat'), which splits the two beams in the vertical plane and then recombines them 
on a common axis. The beams pass through momentum defining collimators and a series 
of four quadrupoles designed to focus the beams at the detector. Finally the two beams 
are again split in the vertical plane and recombined in a second achromat. The layout of 
the beams and detectors is shown schematically in Fig. [TJ 

The beams then enter the fiducial decay volume housed in a 114 m long cylindrical 
vacuum tank with a diameter of 1.92 m, extended to 2.4 m after the first 65 m. Both 
beams follow the same path in the decay volume: their axes coincide within 1 mm, while 
the transverse size of the beams is about 1 cm. With 7 x 10^^ protons incident on the 
target per SPS spill of 4.8 s duration, the positive (negative) beam flux at the entrance 
of the decay volume is 3.8 x 10^ (2.6 x 10"^) particles per pulse, of which 5.7% (4.9%) are 
{K~). The K^/K~ flux ratio is 1.79. The fraction of beam kaons decaying in the 
decay volume at nominal momentum is 22%. 

A detailed description of the NA48 detector can be found in [TU]. The decay volume is 
followed by a magnetic spectrometer housed in a tank filled with helium at nearly atmo- 
spheric pressure, separated from the vacuum tank by a thin (0.31%Xo) Kevlar composite 
window. A thin-walled aluminium beam pipe of 16 cm outer diameter traversing the centre 
of the spectrometer (and all the following detectors) allows the undecayed beam particles 
and the muon halo from decays of beam pions to continue their path in vacuum. The 
spectrometer consists of four drift chambers (DCH): DCHl, DCH2 located upstream, and 
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Figure 1: Schematic lateral view of the NA48/2 beam line (TAXI 7, 18: motorized beam 
dump/collimators used to select the momentum of the and K~ beams; FDFD/DFDF: 
focusing set of quadrupoles, KABESl-3: beam spectrometer stations not used in this 
analysis), decay volume, and detector (DCHl-4: drift chambers, HOD: hodoscope, LKr: 
EM calorimeter, HAG: hadron calorimeter, MUV: muon veto). Note that the vertical 
scales are different in the two parts of the figure. 



DCH3, DCH4 downstream of a dipole magnet. The magnet provides a horizontal trans- 
verse momentum kick Ap = 120 MeV/c for charged particles. Each DCH is composed of 
eight planes of sense wires. The spatial resolution of each DCH is = cfy = 90 /im. The 
nominal spectrometer momentum resolution is (Jp/p = (1.02 © 0.044 ■ p)% {p in GeV/c). 

The magnetic spectrometer is followed by a plastic scintillator hodoscope (HOD) used 
to produce fast trigger signals and to provide precise time measurements of charged parti- 
cles. The hodoscope has a regular octagonal shape with a transverse size of about 2.4 m, 
and consists of a plane of horizontal and a plane of vertical strip-shaped counters. 

The HOD is followed by a liquid krypton electromagnetic calorimeter (LKr) used for 
photon detection and particle identification. It is an almost homogeneous ionization cham- 
ber with an active volume of 7 m^ of liquid krypton, segmented transversally into 13248 
projective cells, 2x2 cm^ each, and with no longitudinal segmentation. The calorime- 
ter is 27Xo deep and has an energy resolution a{E)/E = 0.032//E © 0.09/E © 0.0042 
{E in GeV). Spatial resolution for an isolated electromagnetic shower is (Tx = Cy = 
0A2/\/^ © 0.06 cm {E in GeV) for the transverse coordinates x and y. 

The LKr is followed by a hadronic calorimeter (HAG) and a muon detector (MUV), 
both not used in the present analysis. 



2 Data analysis 

The —>■ TT^e+e" rate is measured relative to the more abundant —>■ tt^tt^, nor- 
malisation channel (where vr^, — > e~^e~'-f is the so called Dalitz decay). The final states 
of the signal and normalisation channels contain identical sets of charged particles. Thus 
electron and pion identification efficiencies, potentially representing a significant source 
of systematic uncertainties, cancel in the first order. 
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Monte Carlo simulation 

In order to compute acceptances for signal, normalisation and background channels, a 
detailed GEANT3-based [H] Monte Carlo (MC) simulation is employed, which includes 
full detector geometry and material description, stray magnetic fields, DCH local ineffi- 
ciencies and misalignment, detailed simulation of the kaon beam line, and time variations 
of the above throughout the running period. 

Event selection 

Three-track vertices (compatible with the topology of K"^ — *■ vr^e^e^ and K"^ — >■ tt^tt^) 
decays) are reconstructed by extrapolation of track segments from the upstream part of 
the spectrometer back into the decay volume, taking into account the measured Earth's 
magnetic field, stray fields due to magnetization of the vacuum tank, and multiple scat- 
tering. 

A large part of the selection is common to the signal and normalisation modes. It 
requires the presence of a vertex satisfying the following criteria. 

• Vertex longitudinal position is inside fiducial decay volume: 2'vcrtcx > -^finai collimator- 

• The vertex tracks are required to be consistent in time (within a 10 ns time window) 
and consistent with the trigger time, to be in DCH, HOD and LKr geometric accep- 
tance, and to have momenta in the range 5 GeV/c < p < 50 GeV/c. Track separa- 
tions are required to exceed 2 cm in the DCHl plane to suppress photon conversions, 
and to exceed 15 cm in the LKr front plane to minimize particle misidentification 
due to shower overlaps. 

• Total charge of the three tracks: Q = ±1. 

• Particle identification is performed using the ratio E/p oi energy deposition of the 
particle in the LKr calorimeter to its momentum measured by the spectrometer. 
The vertex is required to be composed of one vr candidate {E/p < 0.85), and a 
pair of oppositely charged candidates {E/p > 0.95). No discrimination of pions 
against muons is performed. 

If several vertices satisfy the above conditions, the one with the best vertex fit quality is 
considered. The — > vr^e'^'e" candidates are selected by applying the following criteria 
to the reconstructed kinematic variables. 

• 7r^e~*'e~ momentum within the beam nominal range: 54 GeV/c < Ipneel < 66 GeV/c. 

• iT^e'^e~ transverse momentum with respect to the beam trajectory (which is pre- 
cisely measured using the the concurrently acquired Svr^ sample): p"^ < 
0.5 X 10-3 (GeV/c)2. 

• Kinematic suppression of the main background channel vr^vr^ (and other 
minor backgrounds induced by vr^ and vr^^ 4e^ decays) by requiring the e~^e~ 
mass to be above the vr" mass: z = {Mee/ MkY > 0.08, which approximately 
corresponds to Mee > 140 MeV/c^, and leads to an unavoidable loss of ~ 30% of 
the signal sample. 
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Figure 2: (a) Reconstructed spectrum of TT^e~^e~ invariant mass: data (dots) and MC 
simulation (filled area). Note the description of the radiative mass tail by the PHOTOS 
simulation, (b) Reconstructed spectrum of tt^tt^) invariant mass: data (dots) and MC 
simulations of K"^ vr^'^Tr^, signal and K"^ — >■ tt^ii'^u background (filled areas). Signal 
regions are indicated with arrows. 

• 7r^e"'"e~ invariant mass: 470 MeV/c^ < M^ee < 505 MeV/c^. The signal region, 
asymmetric with respect to the nominal kaon mass (Fig. [2^), includes a part of 
the radiative mass tail. The lower mass limit of 470 MeV/c^ corresponds to an 
upper cutoff < 23.1 MeV for the energy of a single directly undetectable soft 
radiative photon, and is chosen as a trade-off between making an inclusive analysis, 
and keeping the background reasonably small. 

No restrictions are applied to the —>■ 7c^e'^e~ sample on the additional energy de- 
position in the LKr calorimeter, the reason being two-fold: 1) to avoid bias due to the 
radiative 7r^e"'"e~7 decays (with predominantly soft photons, which makes the LKr 

response difficult to simulate); 2) to decrease sensitivity to accidental activity. 

For the — > tt^tt^, normalisation mode candidates, the presence of a LKr energy 
deposition cluster (photon candidate) satisfying the following principal criteria is required. 

• Reconstructed cluster energy E > 3 GeV, cluster time consistent with the vertex 
time, sufficient transverse separations from track impact points at the LKr plane 
(i?7r7 > 30 cm, i?e7 > 10 cm). 

• e~^e~7 invariant mass compatible with a vr^j decay: |Mee^ — M^o| < 10 MeV/c^. 

• The same conditions on reconstructed '7r^e^e~7 total and transverse momenta as 
used for 7r^e"'"e~ momentum in the vr^e^e" selection. 

• 7r=^e+e"7 invariant mass: 475 MeV/c^ < M^^ee-y < 510 MeV/c^. 
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Signal sample 

The reconstructed n^e'^e~ invariant mass spectrum is presented in Fig. [2^. The measured 
M^ee resolution is cr,ree = 4.2 MeV/c^, in agreement with MC simulation. The e~^e~ in- 
variant mass resolution computed by MC simulation is = 2.3 MeV/c^. The number of 
K"^ —>■ 7r^e"'"e~(7) decay candidates in the signal region is A^^ee = 7253, out of which 4613 
(2640) are {K~) candidates. The residual background sources after the kinematical 
suppression of the vr^ and tTjjjj decays are the following. 

• Kaon decays with a tc^ in the final state followed by a vr^ — > e^e~7 decay and particle 
misidentification (e^ identified as vr^ or vice versa). Two dominant contributions 
identified with MC simulations are (1) — > tt^tt^, with misidentified and tt^; 
(2) — > TT^e^z/ with a misidentified from the vr^ decay. The two contributions 
are of similar size, which means that about 75% of electron candidates in this sample 
are genuine electrons, and the rest are misidentified pions. 

• Kaon decays with two e"^e~ pairs in the final state resulting from either 7r^(2)) decays 
or external 7 conversions, only a single lepton from each e~^e~ pair constituting the 
best vertex, without particle misidentification involved. All electron candidates in 
this sample are genuine electrons. 

Backgrounds can be reliably estimated and subtracted using the data sample itself. For 
the first of the above background types, the expected kinematic distribution of the lepton 
number violating "same lepton sign" vr'^e^e^ candidates is identical to that of back- 
ground events (up to a negligible acceptance correction). For the second background 
type, the expected sum of distributions of the same lepton sign candidates and the 
TT^e^e^ candidates with the unphysical total charge \Q\ = 3 is similarly identical to 
that of background events. Background contamination is estimated using the sum of the 
numbers of same lepton sign candidates Ngg = 55, and |(5| = 3 candidates Ng^ = 16 as 

{Nss + iV,3)/iV.ee = (l-O ± 0.1stat.)%. 

A cross check of the composition of same lepton sign and \Q\ =3 samples using a 
tighter electron identification (based not only on E/p but also on shower properties and 
track-cluster matching) confirms the expected fractions of genuine electrons among the 
electron candidates. 

Normalisation sample 

The reconstructed 7r^e"'"e~7 invariant mass spectrum is presented in Fig. The number 
of — i> vr^TT^ candidates in the signal region is N2Tr = 1.212 x 10^, which corresponds 
to the number of kaon decays in the fiducial decay volume of about Nk = 1.7 x 10^^. The 
only significant background source is the semileptonic —>■ n^fi^u decay. Its contribu- 
tion is not suppressed by particle identification requirements, since no tt/h separation is 
performed. The background contamination in the signal region is estimated to be 0.15% 
by MC simulation. 

Trigger chain and its efficiency 

Both K"^ TT^e^e^ and ^^^r^ samples are recorded via a two-level trigger chain 

designed to collect the K"^ Stt^ decays. At the first level (LI), the HOD surface is 
logically subdivided into 16 non-overlapping square regions; a coincidence of hits in the 
two planes of the HOD is required to occur in at least two such regions. The second level 
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(L2) is based on a hardware system computing coordinates of hits from DCH drift times, 
and a farm of asynchronous processors performing fast track reconstruction and running a 
selection algorithm, which basically requires the presence of at least two tracks consistent 
with a common origin in the decay volume (closest distance of approach less than 5 cm). 
LI triggers not satisfying this condition are examined further and accepted nevertheless 
if there is a reconstructed track not kinematically compatible with a vr^vr^ decay of a K"^ 
having momentum of 60 GeV/c directed along the beam axis. 

The NA48/2 analysis strategy for non-rare decay modes involves direct measurement of 
the trigger efficiencies using control data samples of downscaled low bias triggers collected 
simultaneously with the main triggers. However direct measurements are not possible for 
the — > 7r^e+e~ events due to very limited sizes of the corresponding control samples. 
Dedicated simulations of LI and L2 performance were used instead. The sources of the 
trigger inefficiency are the local inefficiencies of the HOD and the DCHs, which were 
mapped using special muon runs and the — > tt^tt+tt" data, and included into the 
simulations. The simulated efficiencies and their kinematic dependencies were compared 
against measurements for the abundant K"^ ti^tt^ and tt^'^tt+tt" decays in 

order to validate the simulations; simulated values of trigger efficiencies agree with the 
measurements to a level of a few lO""^ units. 

The simulated values of LI and L2 inefficiencies for the selected vr^vr^ sample 

are £li = 0.37%, = 0.80%. The values of the integral trigger inefficiencies for the 
K"^ — > 7r'^e~^e~ sample depend on the a priori unknown form factor; the corrections 
are applied differentially in bins of dilepton invariant mass Mge. Indicative values of 
inefficiencies for a realistic linear form factor with a slope 6 = 2.3 are Eli = 0.06%, 
£l2 = 0.42%. The K"^ —>■ tt^'^tt^, sample is affected by a larger trigger inefficiency due to the 
smaller invariant masses Mee, which means that the lepton trajectories are geometrically 
closer. Similarly, the inefficiency is enhanced for TT^e~^e~ events at low Mee- 

Theoretical input 

The decay is supposed to proceed through single virtual photon exchange, resulting in a 
spectrum of the z = [Mee/MKY kinematic variable sensitive to the form factor W{z) [2]: 

where Ve = rrie/MK, Tjr = niT^/MK-, and A(a, 6, c) = a? + IP' + — 2ab — 2ac — 26c. On 
the other hand, the spectrum of the angle O-n^e between vr and in the e~^e~ rest frame 
is dT/d9.„e = Csin^^^e, C = const, and is not sensitive to W{z). The expression for 
two-dimensional decay density can be found for instance in pj. 

The following parameterizations of the form factor W{z) are considered in the present 
analysis. 

1. Linear: W{z) = GpM^^fo^l + Sz) with free normalisation and slope (/o, 5)- Decay 
rate and z spectrum are sensitive to |/o|, not to its sign. 

2. Next-to-leading order ChPT |3]: W{z) = GFMl{a+ + h+z) + W^^{z) with free 
parameters (a^., 6+) and an explicitly calculated pion loop term W^'^{z) given in p]. 

3. Combined framework of ChPT and large- A^c QCD [1]: the form factor is parame- 
terized as W{z) = W(w,P,z) with free parameters (w,/?), see also the Appendix. 
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Figure 3: (a) z spectrum of the selected vr^e^e" candidates. Filled area: back- 

ground estimated with same lepton sign and IQI = 3 events, multiplied by a factor of 5 
for visibility, (b) spectrum of the selected candidates (dots) and its comparison to 
MC assuming vector interaction (filled area). 



4. ChPT parameterization [5] involving meson form factors: W(^z) = W{Ma, Mp, z). 
The resonance masses {Ma, Mp) are treated as free parameters in the present anal- 
ysis. 

The goal of the analysis is measuring the form factor parameters in the framework of each 
of the above models, and the computation of the corresponding branching ratios BRi 2,3,4 
by integration of ([T]) normalised to full decay width |14] . 

The Coulomb correction factor is taken into account following [13] (eq. 20). The 
effects of radiative corrections to the 7r^e"'"e~ process are evaluated using the 

PHOTOS [12] simulation of the — ^ 7r^7* TT^e~^e~ decay chain. The results are 
cross-checked using a generalized computation for a multi-body meson decay [13] (eq. 
19). In addition to affecting the values of the model parameters, radiative corrections are 
crucial for the extrapolation of the branching ratio from the limited M^ee (equivalently, 
E^) signal region to the full kinematic region; according to the above models for radiative 
corrections, about 6% of the total —>■ 7i^e^e~{'~f) decay rate fall outside the signal 
region < 23.1 MeV. 

Fitting procedure 

The observed spectra of the data events in z (in the visible region z > 0.08) and cos6'^e 
variables are presented in Fig. [31 The latter spectrum, which contains no information 
about W{z), is compared to the MC expectation, demonstrating a good agreement. 

The values of dTT^ee/dz in the centre of each i-bin of z, which can be directly compared 
to the theoretical predictions ([T]), are then computed as 

(dT^ee/dz). = ^^^^ . ■ ■ - ■ BRiK^ vr^vr") • BR(vr° ). (2) 

iV2^ Ai{l- Ei) Az Tk 

Here Ni and A^^^ are the numbers of observed —>■ n^e'^e~ candidates and background 
events in the i-th bin, is the number of tt'^tt^ events (with background 
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Figure 4: (a) Geometrical acceptances Ai and (b) trigger inefficiencies Ei for the — > 
7r^e~^e~ sample in bins of z, both obtained by MC simulations. The kinematical depen- 
dence of the trigger inefficiency is mainly due to enhancement of LI inefficiency for event 
topologies with two tracks hitting the same HOD segment. The polynomial superimposed 
on Figure (b) is drawn only to guide the eye. 



subtracted), Ai and Ei are geometrical acceptance and trigger inefficiency in the i-th bin 
for the signal sample (computed by MC simulation, both presented in Fig.H]), A2n = 2.95% 
and E2n = 1.17% are those for K"^ tt^'^tt^, events, Az is the z bin width, and is set to 
0.02. The external inputs are the kaon lifetime tk = (1.2380 ± 0.0021) x 10"^ s, and 
branching ratios of the normalisation decay modes BR(i^^ vr^vr") = (20.68 ± 0.13)%, 
BR(7rO,) = (1.198 ±0.032)% [H]. 

While the simulated values of geometric acceptances and trigger efficiencies are subject 
to relatively large systematic uncertainties (discussed in the following section), their ratios 
Ai/A2TT and (1 — Ei)/{1 — E2n) appearing in ([2]) are affected by these uncertainties to a 
lesser extent due to cancellations of most systematic effects. 

The computed values of dF^ee/c^^ vs z and the results of the fits to the four consid- 
ered models are presented in Fig. The corresponding squared form factors |Vr(z)p 
normalized by a condition |iy(0)| = 1 are presented in Fig. \5jp. 

Branching ratios BR(i^^ 7r^e+e~) in the full kinematic range corresponding to each 
model are computed using the measured value of the model parameters, their statistical 
uncertainties, and correlation matrices. In addition, a model-independent branching ratio 
BRjni in the visible kinematic region z > 0.08 is computed by integration of dTj^ee/dz. 
BRmi differs from each of the mo del- dependent BRs computed in the visible z range by 
less than 0.01 x 10~^. The differences between the model-dependent BRs come from the 
region z < 0.08, as can be seen in Fig. [5^. 

Systematic uncertainties 

Particle identification. Imperfect MC description of electron and pion identification 
inefficiencies fe and /^r can bias the result only due to the momentum dependence of fe 
and due to identical charged particle composition, but differing n and e momentum 
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Figure 5: (a) dTT^ee/dz (background subtracted, corrected for trigger efficiency) and fit 
results according to the four considered models, (b) Squared form factors for 
each model normalized by a condition |W^(0)| = 1. 



spectra of the signal and normalisation modes. Inefficiencies were measured with the data 
to be in the ranges 1.6% < fj, < 1.7% and 1.1% < fe < 1.7%, depending in momentum. 
Systematic uncertainties due to these momentum dependencies not perfectly described 
by MC were conservatively estimated assuming momentum- independent MC fe and /^r. 

Beam line description. Despite the careful simulation of the beamline including 
time variations of its parameters, the residual discrepancies of data and MC beam geome- 
tries and spectra may bias the results. To evaluate the related systematic uncertainties, 
variations of the results with respect to variations of cuts on track momenta, LKr clus- 
ter energies, total and transverse momenta of the 7r^e~^e^(7) final states, track distances 
from beam axis in DCH planes, and distances between particle impact points at the LKr 
calorimeter surface were studied. Stability of the results in bins of decay vertex lon- 
gitudinal coordinate was studied as well. The maximum observed variations of the fit 
parameters are taken conservatively as the systematic uncertainties. 

Background subtraction. For background subtraction in the TT^e~^e~ sample the 
spectra of same lepton sign and \Q\ =3 events were used, as described above. This 
method has a limited statistical precision (with an average of ~ 3 background events in a 
bin of z). The uncertainties of the measured parameters due to background subtraction 
were conservatively taken to be equal to the corrections themselves. 

Trigger efficiency. As discussed earlier, the corrections for trigger inefficiencies were 
evaluated by simulations. In particular, LI and L2 corrections for the measured BR have 
similar magnitudes of a few 10~^. No uncertainty was assigned to the LI correction, due 
to relative simplicity of the trigger condition, and the consequent robustness of the simula- 
tion. On the other hand, the uncertainty of the L2 efficiency correction was conservatively 
taken to be equal to the correction itself. 

Radiative corrections. Uncertainties due to the radiative corrections were evaluated 
by variation of the lower 7r^e"'"e~ invariant mass cut. The results show reasonable stability 
with respect to the variation, as can be expected given a good MC description of the 
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Table 1: Corrections and systematic uncertainties (excluding external sources). 



Parameter 


Particle 


Beam 


Background 


Trigger 


Radiative 


Fitting 






simulation 


subtraction 


meniciency 


corr. 


method 


l/ol 


U.UUl 


O.OOd 


0.001 ± 0.001 


n nno i n nno 

0.002 ± 0.002 


n nn a 

0.004 


n nno 

0.003 


r 



0.01 


0.04 


n n i n n 

—0.05 ± 0.05 


nno 1 nno 

—0.03 ± 0.03 


0.04 


n no 

0.03 




0.001 


0.005 


n nn 1 i n nn t 

—0.001 ± 0.001 


n nno i n nno 

—0.002 ± 0.002 


n nn a 

0.004 


n nn /i 

0.004 


"+ 


n nnn 
U.uuy 


U.Uio 


n noi -U n noi 


n ni K 4- n ni K 
U.Uio ± U.Uio 


n ni /I 
U.Ui4 


U.UiU 


w 


0.001 


0.002 


n nno _i_ n nno 
— 0.002 ± U.002 


n nni _t_ n nni 
— U.UOi ± O.OOi 


n nno 
0.002 


n nn 1 
O.OOi 




n Dfi 
u.uu 


n DQ 


_n no -1- n HQ 


_n nfi -4- n Dfi 
— U.UU HI u.uu 


U.UU 


n OA 

U.U^ 


MjGeV 


0.004 


0.009 


0.009 ± 0.009 


0.007 ± 0.007 


0.009 


0.006 


Mb/GeV 


0.002 


0.003 


0.004 ± 0.004 


0.003 ± 0.003 


0.004 


0.002 


BRi_4 X 10^ 


0.02 


0.02 


-0.02 ± 0.02 


-0.01 ± 0.01 


0.02 


0.02 


BR„,i X 10^ 


0.02 


0.02 


-0.01 ± 0.01 


-0.01 ± 0.01 


0.02 


n/a 



7r^e~^e~ radiative mass tail visible in Fig. [2^. Estimation of systematic uncertainties as 
differences of fit results obtained with radiative corrections calculated using the PHOTOS 
simulation and the formalism given in [13] leads to similar results as those obtained with 
the first method. 

Fitting method. Uncertainties due to the fitting procedure were evaluated by com- 
paring the fit results obtained with a bin width Az = 0.01 with those obtained with the 
standard bin width of Az = 0.02. 

External input. Substantial uncertainties arise from the external input, as BR(7r£,) = 
(1.198 ± 0.032)% is experimentally known with a limited relative precision of 2.7% [14] . 
The only parameter not affected by this external uncertainty is the slope 6 of the linear 
form factor describing the shape of the spectrum. 

The applied corrections and the systematic uncertainties, excluding those from exter- 
nal sources, are summarized in Table 1. Uncertainties due to the external input can be 
found in Table 2. 

3 Results and discussion 

The measured value of the model-independent BRmi(2; > 0.08), as well as the parameters 
of the considered models and the corresponding BRs in the full z range, with their statisti- 
cal, systematic, and external uncertainties are presented in Table 2. The 68% confidence 
level contours for the pairs of parameters corresponding to each model are presented 
in Fig. [6l The corresponding correlation coefficients between the model parameters are 
p{\fo\,5) = -0.962, p(a+, 6+) = -0.913, p(w, /3) = 0.999 and p{Ma, Mp) = 0.998. 

The measurement of the model-independent BRini(2; > 0.08), while not involving ex- 
trapolation into the region z < 0.08 depending on the form factor parameterization, still 
relies on the extrapolation into the region of the radiative tail M^ee < 470 MeV/c^. 

Each of the four considered models provides a reasonable fit to the data (x^ of the fits 
are indicated in Fig. [5^), however the linear form factor model leads to the smallest x^- 
The size of the data sample is insufficient to distinguish between the models. 

Earlier measurements of the decay parameters are summarized in Table 3. The present 
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Table 2: Model-independent BRjni(2; > 0.08), and fit results for the considered models. 



BR^i X 10^ = 


2.28 


± 


0.03stat. 


± 


0.04syst. 


± 


0.06ext. = 


2.28 


± 


0.08 


Model (1) 






















l/ol = 


0.531 


± 


0.012stat. 


± 


0.008syst. 


± 


0.007ext. = 


0.531 


± 


0.016 


6 = 


2.32 


± 


- 15stat. 


± 


0.09syst 






2.32 


± 


0.18 


BRi X 10^ = 


3.05 


± 


0.04stat. 


± 


0.05syst. 


± 


0.08ext. = 


3.05 


± 


0.10 


Model (2) 






















a+ = 


-0.578 


± 


0.012stat. 


± 


0.008syst 


± 


0.007ext. = 


-0.578 


± 


0.016 


K = 


-0.779 


± 


0.053stat. 


± 


0.036syst. 


± 


0.017ext. = 


-0.779 


± 


0.066 


BRa X 10^ = 


3.14 


± 


0.04stat. 


± 


0.05syst. 


± 


0.08ext. = 


3.14 


± 


0.10 


Model (3) 






















w = 


0.057 


± 


0.005stat. 


± 


0.004,yst 


± 


O.OOlext. = 


0.057 


± 


0.007 


P = 


3.45 


± 


0.24stat. 


± 


. 1 7svst 


± 


0.05ext. = 


3.45 


± 


0.30 


BR3 X 10^ = 


3.13 


± 


0.04stat. 


± 


0.05syst. 


± 


0.08ext. = 


3.13 


± 


0.10 


Model (4) 

Ma/GeV/c^) = 






















0.974 


± 


0.030stat. 


± 


0.019,yst. 


± 


0.002ext. = 


0.974 


± 


0.035 


M,/(GeV/c2) = 


0.716 


± 


O.Ollstat. 


± 


0.007syst. 


± 


0.002ext. = 


0.716 


± 


0.014 


BR4 X 10^ = 


3.18 


± 


0.04stat. 


± 


0.05syst. 


± 


0.08ext. = 


3.18 


± 


0.10 



measurement of the linear form factor slope 6 is in agreement with the previous ones 
based on 7r+e+e~ [3, [8] and n^fi'^fi' j^Sj samples, and further confirms the 

contradiction of the data to meson dominance models [TB] which predict lower values of 
the slope parameter. The obtained |/o|, a+ and 6+ are in agreement with the only previous 
measurement [8]. The obtained values of w and (3 are in fair agreement with the fit of the 
data points [8] performed by the authors of [1] (the value of Mp = 0.77549 GeV/c^ [13] is 
used in our fit; using Mp = 0.7711 GeV/c^ as in [1] leads to a central value of /3 = 3.27, 
which improves agreement with [4J). The measured parameter Ma is consistent with the 
nominal mass of the ao(980) resonance, while Mp is a few % lower than the nominal p(770) 
mass [13] • 

The branching ratio in the full kinematic range, which is computed as the average 



Table 3: Summary of earlier measurements of the form factor parameters and BR with 
their statistical, systematic, and mo del- dependence uncertainties. Systematic uncertain- 
ties were not evaluated for some measurements; a single statistical uncertainty is presented 
for these cases. 



Decay 
Refs. 



6 

l/ol 
a+ 

b+ 
w 



K: 



[6] 



1.31 ±0.44 ±0.19 



2.14 ±0.13 ±0.15 
0.533 ±0.012 
-0.587 ±0.010 
-0.655 ±0.044 
0.045 ± 0.003 
2.8 ±0.1 



K+ 



^•^3-0.95 



BRxlO^ 2.7 ±0.5 2.75 ±0.23 ±0.13 2.94 ± 0.05 ± 0.13 ± 0.05 
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Figure 6: 68% confidence level contours for the parameters of the four models. Only 
statistical and background subtraction errors are included. 

between the two extremes corresponding to the models (1) and (4), and includes an 
uncertainty due to extrapolation into the inaccessible region z < 0.08, is 

BR = (3.11 ± 0.04stat. ± O.OSsyst. ± 0.08ext. ± O.OT^odei) X 10-' = (3.11 ± 0.12) x 10"'. 

The uncertainty of the above result is correlated with those of the earlier measurements 
(those due to external input and mo del- dependence uncertainties). A comparison to the 
most precise of these measurements by the BNL E865 collaboration [8] using the same 
external input, and taking into account correlation of external uncertainties between the 
two analyses, shows a 1.6a level of agreement. In conclusion, the obtained BR is in 
agreement with the previous measurements [HIEIE]- 

Measurements of the decay parameters and BRs were performed separately for 
and K~ decays. In particular, the branching fractions in the full kinematic range and 
their statistical uncertainties were measured to be 

BR+ = (2.99 ± 0.05stat.) x 10"^ BR" = (3.13 ± 0.08stat.) x 10"^. 

This allows making a first measurement of the direct CP violating asymmetry of K'^ and 
K~ decay rates in the full kinematic range. Considering only the uncorrelated systematic 
uncertainties (those due to background subtraction) between and K~ samples, 

A{K^J = (BR+ - BR-)/(BR+ + BR~) = (-2.2 ± l.Sstat. ± 0.6syst.) x 10~^ 
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The other systematic uncertainties are correlated between the and K" samples, and 
thus cancel in A(i^'^gg). A conservative limit for the charge asymmetry of |A(i^'^gg)| < 
2.1 X 10~^ at 90% CL can be deduced from the above value. This result is compatible 
with CP conservation, however the achieved precision is far from the Standard Model 
expectation |A(/s:4e)l ~ lO'^ [3] and even the SUSY upper limit |A(ir^,J| ~ IQ-^* [HIIIH] 
for the CP violating asymmetry. 

Summary 

From a sample of 7253 K"^ 7r^e'^e~ decay candidates with 1.0% background con- 
tamination, the branching fraction in the full kinematic range has been measured to be 
(3.11 ± 0.12) X 10~^, in agreement and competitive with previous measurements. The 
shape of the form factor which characterizes the decay has been evaluated in the frame- 
work of four models, giving consistent results with previous measurements. The first 
simultaneous observation of both charge kaon decays into n^e~^e~ allowed to establish an 
upper limit for the CP violating asymmetry of and K~ decay rates of 2.1 x 10~^ at 
90% CL. 

Acknowledgements 

We gratefully acknowledge the staff of the CERN SPS accelerator and the beamline for the 
excellent performance of the beam. We also thank the technical staff of the participating 
laboratories, universities and affiliated computing centres for their efforts in operation of 
the experiment and data processing. We are grateful to Samuel Friot, Gino Isidori and 
Victor Pervushin for valuable discussions. 

Appendix 

The paper [1] that formulates one of the considered models does not contain the full 
definition for the form factor fv{z)'- the function appearing in the Eq. (3.19) for 
/y(z) is defined only below the 27r threshold. The expressions for xi^) both below and 
above the threshold, kindly provided by the authors of Ref . [4] , are presented below. Let 
us define r = (Mj^/m^)^. Then for z < 4/r (below the threshold), 

/ N 4 4 1/4 / 1 \ 

^^") = 9-3^ + 3 b'V I^TfTT j ' 

and for z > A/r (above the threshold). 
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